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Abstract 

The small GTPase RhoA is a major regulator of actin reorganization during the formation of stress fibers; thus identifying 
molecules that regulate Rho activity Is necessary for a complete understanding of the mechanisms that determine cell 
contractility. Here, we have Identified Arhgap28 as a Rho GTPase activating protein (RhoGAP) that switches RhoA to Its 
Inactive form. We generated an Arhgap28-LacZ reporter mouse that revealed gene expression in soft tissues at El 2.5, pre- 
bone structures of the limb at El 5.5, and prominent expression restricted mostly to ribs and limb long bones at El 8.5 days 
of development. Expression of recombinant Arhgap28-V5 In human osteosarcoma SaOS-2 cells caused a reduction in the 
basal level of RhoA activation and disruption of actin stress fibers. Extracellular matrix assembly studies using a 3- 
dimensional cell culture system showed that Arhgap28 was upregulated during Rho-dependent assembly of the ECM. Taken 
together, these observations led to the hypothesis that an Arhgap28 knockout mouse model would show a connective 
tissue phenotype, perhaps affecting bone. Arhgap28-null mice were viable and appeared normal, suggesting that there 
could be compensation from other RhoGAPs. Indeed, we showed that expression of Arhgap6 (a closely related RhoGAP) was 
upregulated In Arhgap28-null bone tissue. An upregulation in RhoA expression was also detected suggesting that Arhgap28 
may be able to additionally regulate Rho signaling at a transcriptional level. MIcroarray analyses revealed that Col2al, 
Col9al, Matn3, and Comp that encode extracellular matrix proteins were downregulated in Arhgap28-null bone. Although 
mutations in these genes cause bone dysplasias no bone phenotype was detected in the Arhgap-28 null mice. Together, 
these data suggest that the regulation of Rho by RhoGAPs, including Arhgap28, during the assembly and development of 
mechanically strong tissues Is complex and may involve multiple RhoGAPs. 
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Introduction 

The actin cytoskeleton is fundamental to a wide range of cellular 
functions including cellular contractility, stifFness sensing, tissue 
formation, cell migration and cell polarity but the molecular 
mechanisms are complex and not fuUy understood. Members of 
the family of Rho guanosine triphosphatases (GTPases) are major 
regulators of the assembly of actin-based stress fibers along with 
mammalian diaphanous 1 (mDia) and Rho-associated kinase 
(UOCK) [1-3]. ROCK positively drives the assembly of 
contractile actin stress fibers by directly phosphorylating the 
myosin light chain (MLC), and also by inactivating MLC 
phosphatase [2,4] . Dynamic reorganization of the actin cytoskel- 
eton into stress fibers is essential for fibronectin assembly and is 
regulated by signaling from Rlio GTPases [5,6]. Actomyosin 
contractility is required for the translocation of fibronectin-bound 
integrins in specialized cell-matrix adhesions along actin stress 
fibers, a process that is believed to stretch folded fibronectin dimers 
to facihtate their assembly [7-9] . The mechanism for stress fiber- 
mediated ECM assembly is mechano-sensitive (via cell-matrix 



adhesions) and tightly regulated. For examples, disruption to actin 
polymerization or loss of tension caused the misalignment of 
collagen fibrils in a newly synthesized tendon ECM [10,1 1], and 
targeted ROCK overexpression in the epidermis led to increased 
coUagen deposition and ECM stiffness [12]. While it is well known 
that Rho GTPases regulate actin stress fiber assembly, how they 
are regulated during tissue morphogenesis is less well understood. 

Although Rho is a GTPase its rate of GTP hydrolysis slow. 
Efficient hydrolysis of GTP requires a Rho GTPase activating 
protein (RhoGAP), which accelerates the hydrolytic activity up to 
105-fold [13]. There are over 70 genes encoding proteins that 
contain a RhoGAP domain [14]. This multitude of RhoGAPs is 
thought to ensure signaling specificity, for example, via tissue- 
specific expression, specificity for a single GTPase or signaling 
pathway, or that some RhoGAPs act as scaffold proteins or 
effectors for crosstalk between Rho GTPases and other signaling 
pathways (reviewed by [14]). 

We show here that Arhgap28 is differentially regulated during 
mouse embryonic development. The functions of Arhgap28 have 
not been reported but its differential expression has been listed in a 
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variety of cDNA microarray studies, as summarized in Table S 1 . 
Based on amino acid sequence similarities Arhgap28 is closely 
related to ArhgapG, Arhgap 1 1 a, Arhgap lib, Arhgap 1 8, Arh- 
gap40, DLCl (Arhgap?), DLC2 (Arhgap27) and DLC3 (Arli- 
gap38). Some of these RhoGAPs have been shown to regulate 
actin reorganization. Knockout of DLCl is embryonic lethal at 
El 0.5 of mouse development and examination of fibroblasts 
isolated from E9.5 mouse embryos revealed disrupted stress fibers 
and focal adhesions [15]. Mice lacking functional ArhgapB protein 
are phenotypicaUy normal, despite the fact that ArhgapG is a 
RhoGAP for RhoA and causes the loss of actin stress fibers in 
cultured cells [16]. Arhgap 18 also has specificity for RhoA and 
disrupts actin stress fibers, where knockdown of Arhgap 18 can 
enhance stress fiber formation [17]. On the basis of these studies, 
we hypothesized that Arhgap28 regulates actin stress fiber 
assembly. 

Results 

Sequence alignment predicts that Arligap28 has a 
RhoGAP function similar to Arhgap6 and ArhgapIS 

RhoGAP function is mediated via the RhoGAP domain 
(represented schematically in Figure lA), which enhances hydro- 
lysis of GTP by the target Rho GTPase. Alignment of the 
RhoGAP domains of murine ArhgapG, ArhgapIS and Arhgap28 
showed high homology and confirms the presence of a putative 
catalytic arginine residue (R425) in Arhgap28 (Figure IB), 
suggesting that Arhgap28 has RhoGAP function. No other 
conserved domains were identified. 

Expression of Arhgap28-V5 inhibits RhoA activation and 
actin stress fiber assembly 

To examine if Arhgap28 can regulate RhoA signaling and actin 
stress fiber formation, we created a V5-tagged Arhgap28 
expression clone (Arhgap28-V5; see Figure SI for details) and 
transfected SaOS-2 cells. SaOS-2 cells were selected because they 
form prominent stress fibers when cultured on plastic or glass. 
Expression of Arhgap28-V5 was confirmed in transiently trans- 
fected cells by western blotting using an anti-V5 antibody 
(Figure 2A). Active RhoA was examined by a Rhotekin-GST 
pull-down assay to precipitate GTP-bound Rho. Expression of 
Arhgap28-V5 caused a reduction in active RhoA compared to 
cells treated with the transfection reagent only or when transfected 
with the empty vector (data not shown). In further experiments we 
performed quantitative ELISA assays for active RhoA, Racl and 
Cdc42 in the presence and absence of expressed Arhgap28-V5 
and a mutant Arhgap28R425A-V5 (R425A-V5) with control 
samples transfected with empty vector only. Note that the 
R425A-V5 construct lacks the putative catalytic arginine residue 
in the GAP domain (see Figure 1). We consistently observed 
reduction in the amounts of active RhoA (but not Racl and 
Cdc42) in the presence of Arhgap28-V5. The amounts of active 
RhoA remaining ranged from G2-78% of basal levels of active 
RhoA observed in controls. The results of the experiment that 
showed reduction to 78% is shown in Figure 2B; the expression of 
Arhgap28-V5 but not the R425A-V5 resulted in a significant 
reduction in basal levels of active RhoA with no significant effect 
on basal Racl and Cdc42 activity. The variability in the results 
between experiments was most probably because of differences in 
transfection efficiency between experiments. In separate experi- 
ments we prepared stably transfected SaOS-2 cells, expressing 
Arhgap28-V5. However, for reasons that were unclear to us, the 
expression levels were very low. Overall, the experiments showed 
that Arhgap28-V5 was effective at reducing the levels of RhoA. 



In further experiments, actin stress fibers were examined by 
staining with phaUoidin and immunofluorescence using the anti- 
V5 antibody; cells without staining for V5 exhibited prominent 

stress fibers (Figure 2C). In all cells expressing Arhgap28-V5, as 
identified by staining with V5 antibodies, disrupted actin stress 
fibers were observed. Multiple actin microspikes and membrane 
rufiles were also observed on the edge of these cells (Figure 2C, 
arrows). To test if these morphological changes were attributed to 
the GTP hydrolysis activity mediated by Arhgap28, cells were 
transfected with R425A-V5. All cells that stained positively for V5 
contained [jromincnt stress fibers (Figure 2D). These expression 
studies indicated that Arhgap28 is a RhoGAP for RhoA and its 
expression negatively regulates stress fibers. 

Arhgap28 has a restricted expression pattern through 
embryonic development 

We studied the spatial and temporal activation of Arhgap28 
during embryonic development using an Arhgap28 gene trap 
{Arhgap28^) mouse. The gene trap cassette does not disrupt the 
expression of the wild type Arhgap28 transcript and so the animals 
develop normally (see Figure S2 for details). In this mouse, the 
endogenous promoter of Arhgap28 drives the expression of |3- 
galactosidase, which can be localized by X-gal staining. We 
examined the expression of p-galactosidase in Ell. 5 to E18.5 
embryos. At all time points examined, |3-galactosidase activity had 
a restricted spatial pattern (Figure 3 A). At El 1.5 and E12.5, 
staining was localized to the dorsal region in what appears to 
contain somatic cells (Figure 3 A arrows in lower panels). At El 3.5 
to El 5.5 the staining in the body of the embryos spread to the 
limbs and regions where the ribs are formed. At E18.5, staining 
could be clearly seen in hmb bones and in the dorsal portions of 
the ribs. El 8.5 embryos were sectioned and counter stained with 
Alizarin red. The results showed that cells staining positively for |3- 
galactosidase activity were localized to the calcified portions of 
long bones and ribs (Figure 3B). These data suggested to us that 
Arhgap28 has a role in regulating Rho in the initial stages (from 
E12.5) of ECM assembly throughout the embryo and becomes 
restricted to boney tissues at late stages of embryonic development, 
in the mouse. 

Arhgap28 and related RhoGAPs are upregulated during 
Rho-mediated tissue stiffening in vitro 

To further explore the possibility that Arhgap28 is involved in 
the assembly of ECM, we examined the expression of Arhgap28 in 
a 3D cell culture model of tissue assembly in which embryonic 
fibroblasts deposit and tension a collagen fibril-rich ECM [1 1] . In 
this system, fibroblasts are moved from conventional 2D culture to 
medium containing fibrinogen and thrombin. The formation of a 
fibrin gel occurs within 5 minutes to produce a loose gel in which 
the fibroblasts find themselves suspended. During the next ~10 
days in culture, the fibroblasts replace the fibrin with a collagen/ 
fibronectin-rich matrix [18] and subsequently tension this matrix 
using non-muscle myosin Il-derived forces [19]. Quantitative PCR 
showed that the expression ot Arhgap28 was low in cells cultured 
on plastic (Figure 4A). The expression increased when the cells 
were in fibrin gels but was notably upregulated (15-fold) once the 
tissue constructs had tensioned (day 13; j!><0.05). The expression 
levels continued to increase and were 35-fold higher (compared to 
2D culture) after 14 days of 3D culture under tension (25 days in 
total; j()<0.001; Figure 4A). Arhgap6 and Arhgapl8 were also 
upregulated during tissue construct formation but not to the same 
extent as Arhgap28 (Figure 4A). 



PLOS ONE I www.plosone.org 



2 



September 2014 | Volume 9 | Issue 9 | el 07036 



Arhgap28 Regulates RhoA and Stress Fiber Assembly 




B 

351 

Arhgap6 STSPNTPEPA 
ArhgaplS RIGDLAPQDM 
Arhgap28 ETGDLSVEDM 
Consensus .tgdl.p#din 



PRARRRGAMS VDSITDLDDN 
KKVCSLSLIE LTALYDVLGL 
KKIRHLSLIE LTAFFDAFGI 
kk.r.lslie lta..D..g. 



: of RhoGAP domain4' 

400 

QSRLLEALQL SLPAEA. .QS 
EFKQQKAVKI KTRDSGLFGI 
QLKRNKTERV RGRDNGIFGV 
#. k. .ka rd. g. f g . 



401 -[ AO ] 

Arhgap6 KKEKARDKKL SLNPIYRQ' 
ArhgaplS PLTILLEQDQ RKVPGTR. 
Arhgap28 PLTVLLDNDR KKDPAVK. 
Consensus plt.ll#.d. .k.P..r. 



I 



•[ A 
RLVDSCCQHL EKHG| 
LIFQKLISRI EEGS, 
LVLQKFFQKV EESGI 
l..#k..q.. Ee.gl 




[ 450 
'GSSKKRV 
IPGAAMRI 
SGCTAKV 
g...r! 



451 
Arhgap6 
ArhgaplS 
Arhgap2S 
Consensus 



Al 

RQLREl 
KNLCQ 
KQYRE 



k#lr#|l#a 



r 



] — 

IFDR. 
EAK 
DAR 



-[ 



.GVDVCLEEE HSVHD' 
FYEGTFNWES VKQHD, 
FNADKFKWDK MCHRE, 
f..d.f.W#. ...h# 




500 



501 
Arhgap6 
ArhgaplS 
Arhgap2S 
Consensus 




550 



551 
Arhgap6 
ArhgaplS 
Arhgap2S 
Consensus 




601 
Arhgap6 
ArhgaplS 
Arhgap2S 



]- 



lAWQKMIEN YE, 
MHL. . .LIRY QKIl 
IRL. . .MLKY QKI 
Consensus i.l...$i.y qki| 




HKQK SSD 

IFMCHT LGl| 
FFSRS 



[ 

IEYSVQS sa: 
SSEQRE FE: 
HSDCEE LL: 
. s . . .e 




600 



"I'End of RhoGAP domain 



650 



|PD LQNEl 
TliKF IVIqI 
IsF LITqI 



jISLL ETDPDWDYL LRRKASQSSS 
IKQNI ENQ. .KKERK AMKKLLKKMA 
IRMNE ATMLLKKQLP SMKKLLRRKT 



.f l..#ir..n. et...kk#.. .mkKll 



Figure 1 . Sequence alignment of the RhoGAP domain of Arhgap28-reiated RhoGAPs. A. Schematic of the helical segments and loops in 
the RhoGAP domain. The catalytic arginine residue is highlighted in pink. B. The amino acid sequences of Arhgap6, ArhgaplS and Arhgap28 were 
aligned and the sequences of RhoGAP domain are shown. The helical segments (AO, A, Al, B, C, D, E, F and G) are shown in yellow. Residues 
conserved in all three RhoGAPs are highlighted in red. The catalytic arginine residue is highlighted in pink. 
doi:10.1371/journal.pone.0107036.g001 



To test if Rho signaling is involved in 3D tissue construct 
contraction, fuUy formed tissue constructs (10 days of culture) were 
treated with lysophosphatidic acid (LPA) to activate Rho and the 
formation of stress fibers [20,21] or the ROCK inhibitor, Y27632. 
As shown in Figure 4B, actin polymerization was dramatically 
affected by these treatments; LPA induced thicker, more 
prominent actin stress fibers whereas Y27632 resulted in stress 
fiber shortening and disassembly. 



The 3D tissue constructs recoil when unpinned, which is a 
process driven by non-muscle myosin Il-dependent actomyosin 
contraction [22] . We made use of this feature to develop an assay 
of Rho-myosin Il-dependent contraction of the tissue constructs. 
Thus, constructs were incubated with either LPA or Y27632 and 
the length of the unpinned construct was measured during 30 
minutes. Tissue constructs treated with LPA significandy con- 
tracted to 53% ±0.02 of the original length in 30 minutes. 
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Figure 2. Arhgap28-V5 inhibits RKioA activation and stress fiber formation in SaOS-2 cells. SaOS-2 cells were transfected with empty 
vector or Arhgap28-V5. A. The expression of Arhgap28-V5 was confirmed by western blotting using an antibody to VS. B. Effect of Arhgap28-VS 
expression on the basal activity of RhoA (n = S), Raci (n = 3) and Cdc42 (n = 3). Bars show SEIVl. * indicates significant difference found, p<O.OS. C-D. 
F-actin in cells expressing Arhgap28-VS (C) and Arhgap28R425A-VS (D) was visualized by fluorescence microscopy using anti-VS antibodies and Atto 
488-conjugated phalloidin (representative images from 3 independent transfections). Arrows point to membrane ruffling and F-actin protrusions. 
Bars =2S |im. 

doi:10.1371/journal.pone.0107036.g002 



compared to .59.6% ±0.04 in DMSO-treated tissue constructs ip< 
0.001; Figure 4C). Y27632 treatment significantly inliibited tlie 
contraction of the construct, which contracted to only 71.2%±0.6 
of the original length after 30 minutes (/)<0.001; Figure 4C). The 
Ti/2 for recoil was: 444 s±16 for DMSO-treated control tissue 
constructs, 323 s±7 for LPA treated constructs, and 670 s±3 for 
Y27632 treated constructs (Figure 4D). In addition, expression of 
Arhgap28 and ArhgapS was dramatically downregulated at 6 and 
24 hours after tissue constructs were unpinned (Figure 4E). 
Together, these data show that Arhgap28 is upregulated and 
Rho signaling is active during contraction of newly synthesized 
ECM, and the expression oi Arhgap28 is inversely related to the 
stiffness of the mechanical environment. 

Normal embryonic development of Arhgap28-null mice 

To test further the possible role of Arhgap28 in tissue formation, 
we generated a functional Arhgap28-null mouse by crossing 
Arhgap28^ mice with Cre transgenic mice. LoxP sites flank exons 



7, 8 and 9 of the Arhgap28^' gene. Therefore, Cre recombinase 
would be predicted to produce mice harboring an Arhgap28 del7- 
9 {Arhgap28^"') aUele (Figure S3A). Using PGR readouts with 
genomic DNA we were able to identify wild type, heterozygous 
and null mice (Figure S3B). The absence of exons 7-9 causes a 
frame shift and a smaller transcript that can be detected by RT- 
PCR using primers that span exons 6 to 11 (Figure S3C). Thus, 
transcripts from the defective allele would produce a truncated 
protein lacking the RhoGAP domain. The results showed that 
heterozygous Arhgap28'^"' mice were viable and fertile. Litters 
produced by cross breeding heterozygotes had a close to normal 
Mendelian distribution of genotypes: 30%±9 wild type; 49%±5 
heterozygotes; and 22% ±8 homozygotes (±SEM; from a total of 
37 pups, in 8 litters). Arhgap28'^"' embryos (E15.5) were examined 
histologically and no apparent differences were observed com- 
pared to wild type mice (Figure S3D). 
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Figure 3. Spatial and temporal expression of Arhgap28\n Arhgap2S^' reporter mice. Embryos from time-mated heterozygous Arhgap28^' 
mice were stained with X-gal and genotyped. A. Embryos homozygous for the Arhgap28^' allele were imaged. Wild type embryos acted as controls 
(not shown). Bars =2 mm. B. Homozygous Arhgap28^' embryos from El 8.5 litters were processed for wax embedding and sagittal sections were 
counterstained with Alizarin red S to stain for calcified matrix. Areas of where X-gal-positive cells are shown. Cells stained positive with X-gal are blue 
green and red shows tissues stained with Alizarin red S. Bars =200 |j,m. 
doi:1 0.1 371 /journal.pone.01 07036.g003 



Arhgapd and RhoA are upregulated in Arhgap28-null 
bone tissue 

To confirm the loss of Arhgap28 in bone tissues, RNA was 
isolated from the tibia and fibula of PO wild type and Arhgap28''"' 
mice, and overlapping PCRs were performed (Figure 5A). Full- 
length Arhgap28 was detected in wild type bone tissues and 
overlapping PCRs confirmed that exons 7-9 were absent in the 
Arhgap28 transcript oi Arhgap28^"' bone tissues (Figure 5B). The 
products of these reactions were sequenced to confirm that the 
transcript expressed in Arhgap28^"' bone encoded a 367 amino 
acid-long protein lacking a functional RhoGAP domain (Figure 
S4). 

Due to the lack of an obvious phenotype in heterozygous or null 
mice, we investigated the possibility that other RhoGAPs (e.g. 
ArhgapG and ArhgaplS) were compensating for the absence of 
Arhgap28. Quantitative PGR showed that the truncated Arh- 
gap28'^''' transcript was detectable at low levels in tissues from 
Arhgap28^'^^ mice, which suggests that this truncated transcript 
(lacking the sequences encoding the RhoGAP domain) is at least 
partially stable. However, it was clear that the mutant transcript 
was at significantly lower levels compared to the expression of the 
Arhgap28 transcript in wild type mice (^<0.001; Figure 6A) 
presumably due to nonsense mediated mRNA decay. We also 
detected a significant 2-fold upregulation oiArhgap6 expression in 
Arhgap28^'^ tissues (/)<0.001). Noteworthy, no significant differ- 
ences were found in the expression oi Arhgapl8 in Arhgap28''"' 
tissues (Figure 6A). The data show that in tissues devoid of 
fiinctional Arhgap28, there is compensatory upregulation of 
Arhgap6. We have previously shown that the two RhoGAPs were 
upregulated during ECM assembly and tensioning and were 
downregulated in response to lack of tension, which indicate that 
the Arhgap6-Arhgap28 pair might co-regulate the same Rho 
signaling pathway for actin reorganization. We also examined the 



expression of Rho GTPases. A small but significant upregulation 
in RhoA expression, but not in expression of Racl, Cdc42 or 
RhoQj was found in Arhgap28^''' bone tissues compared to wild 
type (p<O.Ob; Figure 6B). Again, these data were indicative of a 
role of Arhgap28 in RhoA regulation. 

IVlicroarray comparison of gene expression between 
bone tissues of wild type and ArhgaplS''^' mice 

While Arhgap6 expression was upregulated in Arhgap28^''' 
bone tissues there was also a significant increase in RhoA gene 
expression, which suggests that Rho signaling might still be altered 
in Arhgap28'^'^ tissues despite ArhgapS compensation. To 
determine if global changes in gene expression occurred, 
microarray analyses were performed in which gene expression 
from bone tissues (tibia and fibula) of wild type and Arhgap28^''' 
mice were compared. The integrity of the RNA samples and the 
microarray readouts were analyzed (see Figure S5 for details). A 
total of 45037 probe sets were detected of which 363 of these 
showed significant differential expression of & ± 2-fold difference 
(^<0.05). A heat map generated by hierarchical clustering of all 
the probe sets, based on similarities in the expression level and the 
expression profile is shown in Figure S5D. 

As a cautionary note, three probe sets detected Arhgap28^''^ 
expression in bone from knockout animals (see Table 1). This was 
not unexpected because a 3 ' microarray was used and the mutant 
Arhgap28^^ transcript contains the endogenous 3 '-end of the 
Arhgap28 transcript (see Figure 5). In contrast to the qPCR 
analyses performed previously, readouts for ArhgapS did not show 
a fold change of ^±2 in Arhgap28^''' bone (Table 1). Expression 
of ArhgaplS and Rho GTPases (RhoA, Racl, Cdc42 and RhoQ) 
also showed no differential expression in Arhgap28^'' bone. The 
discrepancy between the microarray and the qPCR analyses could 
be because qPCR is a more quantitative and sensitive method of 
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Figure 4. Arhgap28\s upregulated during ECM assembly and downregulated during Rho-dependent tissue recoil. Primary embryonic 
chick fibroblasts were seeded into a fibrin gel for the formation of 3D cell cultures containing aligned ECM. A. The expression of Arhgap6, ArhgapIS, 
Arhgap28, Collal, and Fn were quantified by qPCR (n = 3). B. 3D tissue constructs were treated with 5 |jM LPA or 10 |iM Y27632 for 2 hours in serum- 
free conditions and changes to the actin cytoskeleton were observed by phalloidin staining. Bars =50 |im. Arrows indicate the alignment of the 
tissue. C. The recoil baseline, which is the calculated length of the construct at the end of the exponential fit (see Materials and Methods for details). 
D. The half-time, which is the time required to achieve half of the recoil baseline. E. The expression of ArhgapS, ArhgapIS and Arhgap2S were 
quantified by qPCR during tissue recoil (n = 3). Fold changes in gene expression were normalized to Gapdh [2^^'^'^^ values). Bars show SEM. *** and * 
indicate significant differences found, p<0.001 and p<0.05, respectively, one way ANOVA. 
doi:10.1371/journal.pone.0107036.g004 



detecting gene expression. Interestingly, the actin genes detected 
by this microarray study, Actal and Add, were both significantly 
downregulated (~ 3-fold) in the bone tissue of Arhgap28''"' mice 
{q<Q.Q5; Table 1). The collagen genes expressed in bone tissues, 
Collal, Colla2 and CollOal, were not significantly different in 
Arhgap28-null bones. 

DAVID online tool was used to perform gene ontology analyses 
on probe sets showing a differential expression &± 2-fold. The 
analysis showed that 57 probe sets detected genes that were 
downregulated and 306 probe sets detected genes that were 



upregulated, in the mutant samples. DAVID identified 10 
functional annotation clusters in the genes that were downregu- 
lated. The top 3 annotation clusters and the genes that are over- 
represented in these clusters are listed in Table 2. All three clusters 
contained genes indicative of a cartilage ECM, for example, 
Col2al, Col9al , Haplnl (hyaluronan and proteoglycan link 
protein 1), Matnl (matriUin 1) and Matn3, and some genes that 
are linked to negatively regulating bone mass, Agtr2 (angiotensin 
II receptor, type 2) and Dlhl, which encodes a transmembrane 
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Figure 5. Comparison of tKie fuli-iengtli transcripts of wild type and Arhgap2Sf'"' a\\&\es. A. Overlapping products of PCR reactions 1,4, 6-7, 
9-14 by specific primers (numbers in gray) are sequenced. ATG and TAG indicates the start and stop codons, respectively. Diagram not drawn to 
scale. B. RNA was isolated from the tibia and fibular of PC mice, and the 10 overlapping RT-PCRs were performed and visualized by gel 
electrophoresis. Products were subsequently purified and confirmed by DNA sequencing. Reactions in the green boxes will not produce a product 
with cDNA from Arhgap28^''' (see Figure S4). 
doi:1 0.1 371 /journal.pone.Ol 07036.g005 



protein called delta-like 1 homolog. These data suggest that a bone 
phenotype could still be expected upon closer analysis. 

For genes that were upregulated in Arhgap28-nuU bone tissue, 
65 annotation clusters were identified and the top 3 clusters are 
listed in Table 3. The genes that were over-represented in the first 
cluster were genes involved in targeting proteins for ubiquitination. 
For example, this cluster contained 13 E3 ubiquitin protein ligases, 
including March3 (membrane-associated ring finger (C3HC4) 3), 
MarchS, Mybp2 (MYC-binding protein 2); 3 peptidases; and a 
gene, Psmdl4, which encodes a regulatory subunit of 26S 
proteasome (see Table 3 for full list of genes). The second cluster 
contained genes that promote actin polymerization, including 
Pafahlhl (platelet-activating factor acetyl-hydrolase lb), Tmodl 
(tropomodulin 1) and Diap3 (also known as mDia2)\ actin 
nucleation [Spirel); and genes involved in linking the actin 
cytoskeleton to the plasma membrane including Utrn (utrophin) 



and Spnal (spectrin). The third cluster of over-represented genes 
that were upregulated in Arhgap28'^''' bone tissue contained zinc 
finger proteins (Rhm5, Neil3 and Fus), which bind DNA or RNA. 
To summarize, the GO analyses showed that loss of Arhgap28- 
mediated RhoA signaling causes: (i) down-regulation of cartilage 
ECM genes; (ii) up-regulation in genes involved in targeting 
proteins for degradation by ubiquitination; and (iii) up-regulation 
of genes involved in anchorage of the actin cytoskeleton to the 
plasma membrane. Whether or not Arhgap28 regulation of Rho 
signaling is involved functionally in cartilage homeostasis, protein 
degradation of anchorage of actin to the plasma membrane wiU 
require further investigation. However, the results presented in 
Figure 2 are indicative of a role for Arhgap28 in actin stress fiber 
polymerization at the plasma membrane. 
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Figure 6. Expression of RI10 GAPs and Rho GTPases in 

Arhgap2Sf'"' bone tissues. RNA was isolated from bone tissues (from 
tibia and fibula) of PO wild type and Arhgap28^''' mice. A and B. The 
expression of genes encoding Arhgap28 and related RhoGAPs, Arhgap6 
and Arhgap18 (panel A); and Rho GTPases, RhoA, Raci, Cdc42 and 
RhoQ (panel B) was quantified by qPCR (n = 3). Fold changes in gene 
expression normalized to Gapdh (2"'^*'-' values). Bars show SEM. *** 
and * indicate significant differences found when compared to wild 
type, p<0.001 and p<0.05, respectively. 
doi:10.1371/journal.pone.0107036.g006 

Analyses of bone length in mature mice show no 
evidence of bone dysplasia 

The results of the GO analyses of Arhgap28-null bone tissues 
showing down-regulation of genes encoding cartilage ECM 
molecules and the fact that Arhgap28 expression was localized 
to pre-bone structures in embryos, prompted us to explore the 
possibility of a bone phenotype in the knockout mice, especially 
because mutations in Col2al , Col9al, Matnl and Matn3 are 
associated with bone dysplasia (for reviews see [23-25]). Thus, we 
examined the skeletons of mice by X-ray at 1 0 weeks of age, when 
the skeletal growth has reached maturity. Skeletal analyses 
revealed no significant differences in the lengths of femur and 
tibia between wild type, heterozygotes and Arhgap28'^''' mice 
(Figure 7). As a measure of intramembranous ossification, we 
measured the intercanthal distance and found no significant 
differences between wild type, heterozygous and Arhgap28-null 
animals (Figure 7). Together, these data suggest that although 
genes associated with bone dysplasia were downregulated in the 



Arhgap28' " mouse there is no bone length phenotype, presum- 
ably because of compensation from ArhgapG. 

Discussion 

The mechanical stiffness of musculoskeletal tissues is directly 
related to the organization of coUagen fibrils in the ECM. For 
examples, the strongest tensile tissues such as tendon and ligament 
have coUagen fibrils arranged in parallel register, which presum- 
ably is the best organization to resist uniaxial force; in bone the 
coUagen fibrils provide a template for mineralization. Further- 
more, the collagen fibrils are pre-stressed by cells to ensure that 
tissues can respond directly to applied forces. It is poorly 
understood how the tissue-specific arrangement and pre-stressiiig 
of collagen fibrils is achieved. Tissue stiffening requires cellular 
contraction via actin stress fibers, which are regulated by Rho 
GTPases. How these Rho GTPases are regulated during tissue 
development to ensure that the ECM is at optimal stiffness is 
unknown. 

Here, we show that Arhgap28 is activated in bone tissues before 
birth and during the assembly of a stiff ECM. Experiments 
expressing Arhgap28-V5 suggest that Arhgap28 is a negative 
regulator of RhoA and actin stress fiber formation. Arhgap28- 
deficiency does not appear to affect bone development, which is 
most likely due to functional redundancy between Arhgap28 and a 
closely related RhoGAP, ArhgapG. It will be important in future 
work to investigate the mechanisms of how Arhgap28-regulated 
actin contractility determines stiffening of the ECM and to 
understand how RhoGAPs crosstalk regulates Rho and actin 
remodeling within developing musculoskeletal tissues. 

Cellular tension is generated by actin stress fibers [26]. In this 
context, the closely related ArhgapG and Arhgapl8 regulate the 
formation of actin stress fibers via RhoA [16,17] and in this study, 
expression of Arhgap28-V5 caused similar effects. Sustained RhoA 
activity has inhibitory effects on Racl- and Cdc42-activated 
lamellipodia and fdopodia formation [27,28], which helps explain 
the appearance of actin microspikes and membrane ruffles in 
Arhgap28-overexpressing cells although there was no detectable 
activity of Arhgap28 against Racl and CDC42 in the assays used 
here. 

Cells respond to stiff extracellular matrices via Rho-activated 
actin stress fibers [29-31] and overactive RhoA signaling is linked 
to cancer [see 32 for review] . Here, loss of RlioGAPs, such as 
DLCl, is associated with cancer [33], which suggest that RhoA 
signaling affects cell fate. Overexpressioii of RhoA and Rho 
GTPases are linked to cancer and the stability oiRhoA transcripts 
in cancer cells has been shown to be a result of altered 
polyadenylation signals [34], further suggesting why it would be 
biologically important to have more than one RhoGAP regulating 
Rho GTPases. 

There are a few studies that show that the differentiation of 
mesenchymal stem cells into an osteogenic lineage can be 
influenced by Rho/ROCK signaling, for examples see [35-37]. 
Surprisingly, mice expressing a dominant-negative RhoA have a 
bone sclerotic phenotype, which suggests that lack of RhoA 
signaling enhances mineralization [38]. These disagreements 
linking Rho/ROCK and bone development could be due to the 
need for a balanced signal and crosstalk between Rho/ROCK and 
growth factors rather an aU-or-nothing response. RhoA expression 
was also significantly upregulated in bone tissues of Arhgap28-null 
mice. RhoA expression is activated by transcription factor Myc 
[39], but Myc was not differentially regulated in the microarray 
comparison between wild type and Arhgap28^''^ bone tissues. 
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Table 1. Expression of genes of interest in Arhgap2ff'^' bone tissues. 





Gene 


Gene symbol wild type 


Arhgapllf"' 


Fold change 


q value 


RhoGAP genes 


Arhgap28 


Arhgap28 


271.7 


110.3 


-2.5 


0.037 






82.5 


76.5 


-1.1 


0.345 






67.1 


64.2 


-1 


0.466 


Arhgap6 


Arhgapd 


62.2 


62.1 


-1 


0.656 






271.4 


275.8 


1 


0.607 






20.2 


24.7 


1.2 


0.257 


ArhgaplS 


ArhgaplS 


488.6 


397.5 


-1.2 


0.106 






82.2 


73.2 


-1.1 


0.356 


Deleted in liver cancer 1 


DLCl 


48.6 


52.3 


1.1 


0.413 


Rho GTPases genes 


RhoA 


RhoA 


5112.5 


5060.3 


-1 


0.627 






4777.4 


4752.9 


-1 


0.643 






12.9 


16.3 


1.3 


0.099 


Raci 


Racl 


3672.1 


3349.2 


-1.1 


0.254 






3591.8 


3611.5 


1 


0.645 






14.5 


18.0 


1.2 


0.085 


Cdc42 


Cdc42 


239.8 


190.0 


-1.3 


0.197 






4229.0 


3745.1 


-1.1 


0.191 






3386.3 


3059.5 


-1.1 


0.352 






8834.9 


8114.0 


-1.1 


0.27 


RhoQ 


RhoQ 


2126.4 


2001.0 


-1.1 


0.52 


actin genes 


alpha-skeletal actin 


Actal 


334.1 


94.0 


-3.6 


0.015 


alpha-cardiac-actin 


Acta 


687.1 


223.8 


-3.1 


0.046 


bone collagen genes 


collagen, type 1, alpha 1 


Collal 


10217.8 


10997.6 


1.1 


0.413 






10473.8 


11427.3 


1.1 


0.528 


collagen, type 1, alpha 2 


Collla2 


321.6 


767.0 


2.4 


0.163 






23205.8 


24471.2 


1.1 


0.422 






19122.7 


1 024.8 


1.1 


0.444 


collagen, type X, alpha 1 


Col) Oat 


1917.6 


1 024.8 


-1.9 


0.056 



Comparison of gene expression in wild type and ArhgaplS bone tissues. Mean intensities from hybridization of triplicate samples to probe set(s) for the genes listed. 
doi:1 0.1 371 /journal.pone.Ol 07036.t001 



Due to the large family of mammalian RhoGAPs, it was 
predicted that Arhgap28-deficiency might be compensated for by 
functional redundancy between Arhgap28 and other RhoGAPs. 
Indeed, absence of ArhgapG in mice does not cause an overt 
phenotype, presumably because of the observed compensatory 
mechanisms [16]. It was surprising in our studies thai ArhgaplS 
was not also upregulated because both Arhgap6 and ArhgaplS 
have been shown to negatively regulate RhoA and actin stress 
fibers. These results reveal the potential of a novel co-regulatory 
mechanism for RhoA signaling and actin stress fibers by Arhgap6 
and Arhgap28. How these RhoGAPs are activated during the 
patterning of ECM is unknown. Matrix assembly and detection of 
the ECM occurs at cell-matrix adhesion sites, therefore it is 
hypothesized that Arhgap28 is activated by signals downstream of 
cell-matrix adhesions in similar mechanisms described for other 
RhoGAPs. For examples, upon activation of integrin pi, pl90- 
RlioGAP is activated via tyrosine phosphorylation at the N- 
terminus by Src [40], and more recently, the identification of 



protein-protein interaction domains within DLCl suggest direc- 
tion interaction between DLC 1 and cell-matrix adhesion proteins 
tensin, talin and focal adhesion kinase [41,42]. 

Importantiy, although there was no overt phenotype in the 
Arhgap28^"' mice (perhaps due to compensation by Arhgap6) the 
microarray study revealed downregulation in Actal and Actcl 
genes, which encode isoforms of actin found most abundandy in 
muscle. Downregulation in both these genes has also been 
observed in a microarray study of chondrocyte differentiation in 
vitro [43] and downregulation oi Actal was observed in bone in 
response to mechanical loading [44] , which suggest of change in 
the cells oi Arhgap28^'' bone tissues. Gene ontology analyses also 
revealed that genes that encode ECM molecules associated with a 
cartilage tissue (including Collal, Col9al , Matn3 and Comp) were 
downregulated in Arhgap28-null bone tissues. Mutations in these 
genes are associated with bone dysplasia [for reviews, see 
23,24,25]. Bone tissues from ArhgaplS'''' mice also showed an 
upregulation of genes involved in ubiquitination and actin 
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Table 2. Top 3 annotation clusters of genes downregulated in Arhgap28 bone. 





Cluster components'^ 


p value^ 


Matched 


genes^ Gene names (if annotated) of corresponding probe IDs in the GO list" 


Annotation cluster 7, Enrichment score 6.72 


SP_PIR_KEYWORDS secreted 


4.13E-11 


21 


Egfl6, Haplnl, Mial, Sfrp4, Ccl6, Col2a1, Lectl, Ptn, MatnS, Ecrg4, Matn3, Haplnl, 
Ccll2, Comp, CtqtnfS, Angptll, Haplnl, Nts, MfapS, Haplnl Epyc, Ecrg4, Mmrnl, 
Cmal, Col9al, CxclU 


GO: 0005576- extracellular region 


6.53E-10 


22 


Egfl6, Haplnl, Mial, Sfrp4, Ccl6, Col2al, Lectl, Ptn, Matn3, Ecrg4, Haplnl, Matn3, 
Ccll2, Comp OqtnfS, Angptll, Haplnl, Nts, MfapS, Matnl, Haplnl, Epyc, Ecrg4, 
Mmrnl, Cmal, Col9al, CxclU 


SP_PIR_KEYWORDS signal 


1 .06E-07 


24 


Egfl6, Mial, Sfrp4, Ccl6, Col2al, Snore, DIkl, Ptn, MatnS, Ecrg4, MatnS, Cell 2, Comp, 
ClqtnfS, CdSOOId, Angptll, Nts, MfapS, Matnl, Cpa3, Epyc, Ecrg4, Mmrnl, Cmal, 
Col9al, CxclU 


GO: 0044421 ~ extracellular region part 


1 .57E-07 


14 


Haplnl, Haplnl, Angptll, Egfl6, DIkl, MfapS, Matnl, Col2al, Haplnl, Epyc, Ptn, Matn3, 
Haplnl, MatnS, Col9al, CclU, CxclU, Comp 


UP_SEQ_FEATURE signal peptide 


7.64E-07 


24 


Egfl6, Mial, Sfrp4, Ccl6, Col2al, Snore, DIkl, Ptn, MatnS, Ecrg4, MatnS, CclU, Comp, 
ClqtnfS, CdSOOId, Angptll, Nts, MfapS, Matnl, CpaS, Epyc Ecrg4, Mmrnl, Cmal, 
Col9a1, CxclU 


UP_SEQ_FEATURE disufide bond 


1.86E-06 


21 


Egfl6, Haplnl, Mial, Sfrp4, Ccl6, Lectl, DIkl, Ptn, MatnS, MatnS, Haplnl, CclU, Comp, 
Agtr2, CdSOOId, Angptll, Haplnl, Matnl, CpaS, Haplnl, Epyc Mmrnl, Cmal, Col9al, 
CxclU 


Annotation cluster 2, Enrictiment score 3.78 


GO: 0005578- proteinaceous EClVl 


1 .48E-07 


10 


Egfl6, Haplnl, Haplnl, MfapS, Matnl, Col2al, Haplnl, Epyc, Ptn, MatnS, Matn4, 
Haplnl, Col9al, Comp 


GO: 0031012- extracellular matrix 


2.07E-07 


10 


Egfl6, Haplnl, Haplnl, MfapS, Matnl, Col2al, Haplnl, Epyc, Ptn, MatnS, Matn4, 
Haplnl, Col9al, Comp 


SP_PIR_KEYWORDS extracellular matrix 


1 .48E-05 


7 


Haplnl, Haplnl, Egfl6, Haplnl, MfapS, Col9al, Comp Col2al, Haplnl, Epyc 


GO: 0007155- cell adhesion 


7.13E-03 


6 


Haplnl, Haplnl, Egfl6, Mial, Haplnl, Col9al, Comp Col2al, Haplnl 


GO: 0022610- biological adhesion 


7.19E-03 


6 


Haplnl, Haplnl, Egfl6, Mial, Haplnl, Col9al, Comp, Col2al, Haplnl 


Annotation cluster 3, Enrichment score 2.24 


SP_P1R_KEYW0RDS egf-like domain 


2.32E-04 


6 


Egfl6, MatnS, Mmrnl, Comp, Matnl, DIkl, MatnS 


Annotation cluster analysis of probe sets detecting a significant fold change greater than —2 from wild type to Arhgap28'^^' produced 1 0 clusters. The top 3 annotation 
clusters with the highest enrichment score are listed here. A. The top gene ontology components of the cluster. B. The statistical significance of this grouping where the 



lower the score the more unlikely this clustering is due to chance. C. The number of probe sets that recognize genes contributing to the GO term. D. List of gene names 
of the Affymetrix Mouse Genome 430 2.0 array probe IDs. 
doi:1 0.1 371 /journal.pone.Ol 07036.t002 



reorganization compared to wild type. Most of the genes enriched 
in tlie top cluster were E3 ubiquitin ligases, which mediate the 
specificity of the ubiquitination pathway [reviewed by 45]. Tlie 
target for many of these E3 ligases is unknown. The genes that 
stood out include Hacel, which catalyzes the ubiquitination of 
Racl [46], Mychp2, which can bind Myc [47] but whether or not 
it is inhibitory or enhancing for RhoA transcription is unknown. 
The differential regulation of these genes is indicative of a response 
caused by Arhgap28-deficiencY, which might be independent of 
RhoA signaling if ArhgapG compensates for Arhgap28. 

In conclusion, we describe experiments that point to Arhgap28 
being a functional negative regulator of RhoA during the 
formation of actin stress fibers in cells of mesenchymal origin. 
However, knockout experiments in mice indicated that Arhgap28 
is not autonomous but is functionally redundant in the presence of 
other RhoGAPs including ArhgapG and possibly ArhgaplS. Such 
redundancy presumably underpins the cell's abihty to generate 
actin stress fibers during tissue formation. 

Materials and Methods 

Ethics 

The care and use of all mice in this study was carried out in 
accordance with UK Home Office regulations, UK Animals 
(Scientific Procedures) Act of 1986 under the UK Home Office 



licence (PPL 40/3485). No experimental procedures were 
performed on live animals. All animals were sacrificed by Schedule 
1 cervical dislocation by trained personnel, and all efforts were 
made to minimize suffering. 

V5-His-tagged nnurine Arhgap28 and Arhgap28M2SA 
expression clones 

For the generation of a C-terminus V5-His-fusion Arhgap28 
clone, the sequence encoding V5-His (GKPIPNPLLGLDST- 
(His)6) was introduced immediately 5 ' of the stop codon of the 
endogenous, larger Arhgap28 isoform cDNA. The cDNA clone of 
the short Mus mmculus Arhgap28 isoform (BC066788.1; Source 
Bioscience Geneservice) was digested with EcoRI and BsmBI 
(New England Biolabs, Hertfordshire, UK) following the manu- 
facturer's instructions to produce a 1123 bp fragment of the 5' 
portion of the Arhgap28 clone. For the preparation of the 3' 
portion of endogenous Arhgap28 transcript, a 1209 bp fragment 
was amplified from mouse fibroblast cDNA (made with random 
hexamers) using a forward primer 5' of a unique endogenous 
BsmBI restriction recognition site 3'-AAG ATT TGG GTT GAG 
CGA GAG G-5' and a reverse primer that created a new BstBI 
site 3'- AAT TGG AAG GGG TTG ATG ACC G-5' immediately 
before the stop codon. iProof High-Fidelity DNA Polymerase (Bio 
Rad Laboratories) was used following the manufacturer's protocol. 
The PGR fragment was sequenced and cut with BsmBI and BstBI 
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Table 3. Top 3 annotation clusters of genes upregulated in Arhgap28 bone. 





Cluster components'^ 


p value" 


Matched genes^ 


Gene names (If annotated) of corresponding probe IDs in the GO list'' 


Annotation cluster 1, Enrichment score 5.02 


SP_PIR_KEYWORDS ubi conjugation 
pathway 


1 .47E- 


-06 


19 


March3, Cbib, Herc4, Mycbp2, Ranbp2, MarchS, Radl8, Fbxl7, Mkml, SpopI, Usp7, 
MarchS, Phrl, Hipl, FbxoSO, Hacel, MarchS, Ube2o, Mycbp2, Dtl, Huwel, Mkrnl, 
PsmdM, Usp25, Herd 


GO: 0030163— protein catabolic process 


1 .64E- 


-06 


21 


MarchS, Cbib, Herc4, Mycbp2, Ranbp2, MarchS, RadlS, Fbxl7, Mkrnl, SpopI, Usp7, 
MarchS, Phrl, Hipl, FbxoSO, UspS2, Hacel, MarchS, Ybey, UspS2, Btrc, Ube2o, 
Mycbp2, Dtl, TtllS, Huwel, Mkrnl, PsmdU, Usp2S, Herd 


GO: 0019941- modification-dependent 
protein catabolic process 


1.72E- 


-06 


20 


MarchS, Cbib, Herc4, Mycbp2, Ranbp2, MarchS, RadW, Fbxl7, Mkrnl, SpopI, Usp7, 
MarchS, Phrl, Hipl, FbxoSO, UspS2, Hacel, MarchS, Ybey UspS2, Ube2o, Mycbp2, 
Dtl, TtllS, Huwel, Mkrnl, PsmdM, Usp2S, Herd 


GO: 0009057- macromolecule 
catabolic process 


5.25E- 


-06 


22 


MarchS, Cbib, Herc4, Mycbp2, Ranbp2, MarchS, RadW, Fbxl7, Mkrnl, SpopI, Usp7, 
MarchS, Phrl, PanS, Hipl, FbxoSO, UspS2, Hacel, MarchS, Ybey UspS2, Btrc, Ube2o, 
Mycbp2, Dtl, TtllS, Huwel, Mkml, PsmdU, Usp2S, Herd 


SP_PIR_KEYWORDS ligase 


6.21E- 


-06 


11 


Hacel, MarchS, Cbib, MarchS, Herc4, Mycbp2, Btrc, MarchS, RadlS, Ube2o, Mkrnl, 
Mycbp2, Huwel, Mkrnl, MarchS, Phrl, Hipl 


GO: 0006508- proteolysis 


2.91E- 


-03 


22 


Cbib, MarchS, Herc4, Mycbp2, Ranbp2, MarchS, RadlS, Fbxl7, Mkrnl, SpopI, Usp7, 
MarchS, Phrl, Hipl, FbxoSO, Usp32, Hacel, MarchS, Ybey UspS2, Ube2o, Mycbp2, 
Dtl, TtllS, Huwel, Mkrnl, PsmdU. Kel, Metap2, Usp2S, Herd 


Annotation cluster 2, Enrictiment score 3.38 


GO: 0003779- actin binding 


1.03E- 


-04 


13 


MyolD, Spnal, MyolD, Fpb4.1, Spnal, Spirel, Spirel, Utrn, Fhdcl, Trpm7, TmodI, 
Mtssl, Epb4.9, Epb4.1, Slc4al, Epb4.1, DiapS, Addl, Fhdcl 


GO: 0030036- actin cytoskeleton 
organization 


3.71E- 


-04 


9 


MyolD, Pafahlbl, Trpm7, Rictor, Fhdcl, Tmodl, Mtssl, Epb4.1, Sorbsl, Epb4.1, 
Sorbsl, DiapS, Fhdcl 


GO: 0007010— cytoskeleton organization 


7.02E- 


-04 


12 


MyolD, Pafahlbl, Stradb Trpm7, Rictor, Fhdcl, Tmodl, Mtssl, Epb4.9, Epb4.1, 
Sorbsl, Epb4.1, Sorbsl, DiapS, Cenpe, Fhdcl 


GO: 0007155- cell adhesion 


8.28E- 


-04 


14 


MyolD, Spnal, MyolD, Epb4.l, Pafahlbl, Spnal, Spirel, Spirel, Utrn, Fhdcl, 
Trpm7, Tmodl, Mtssl, Epb4.9, Epb4.l, Slc4al, Epb4.l, DiapS, Fhdcl, Addl 


Annotation cluster 3, Enrictiment score 3. 15 


UP_SEQ_FEATURE zinc finger region: 
RanBP2-type 


1.09E' 


-02 


3 


RbmS, Neils, Fus 



Annotation cluster analysis of probe sets detecting a significant fold change greater than 2 from wild type to Arhgap28 ^ produced 65 clusters. The top 3 annotation 
clusters with the highest enrichment score are listed here. A. The top gene ontology components of the cluster. B. The statistical significance of this grouping where the 
lower the score the more unlikely this clustering is due to chance. C. The number of probe sets that recognize genes contributing to the GO term. D. List of gene names 
of the Affymetrix Mouse Genome 430 2.0 array probe IDs. 
doi:l 0.1 371/journal.pone.01 07036.t003 



restriction enzymes (New England Biolabs). pcDNA6/V5-His C 
vector was prepared by digestion with EcoRI and BstBI restriction 
enzymes (New England Biolabs). The three cut fi-agments were 
purified using QIAquick Gel Extraction kit (Qiagen) and ligated 
using T4 DNA Ligase (New England Biolabs). Arhgap28R42SA 
clone containing a codon change from CGA to GCA (Genscript) 
was subcloned into pcDNA6/V5-His C between EcoRI and BstBI 
sites. DNA was purified from each culture using QIAprep Spin 
Miniprep kit (Qiagen) following the manufacturer's protocol. 

Transient transfections and GTPase activation assays 

NIH3T3 fibroblasts were cultured in DMEM, supplemented 
with 10% (v/v) FCS, 2 mM L-glutamine, 10 000 U/ml penicillin 
and 10 mg/ml streptomycin. SaOS-2 osteosarcoma cells were 
cultured in RPIM-1630, supplemented with 10% (v/v) FCS, 
2 mM L-glutamine, 10 000 U/nil penicillin and 10 mg/ml 
streptomycin. Cells were cultured at 37°C, 5% CO2. DNA was 
transfected using Lipofectamine 2000 (Invitrogen) according to 
manufacturer's protocol. Cells were analyzed 24 hours after 
transfection by indirect immunofluorescence, active GTPase G- 
LISA assays or western blotting. Expression of V5-tagged 
constructs was confirmed by western blotting using a mouse 
monoclonal anti-V5 epitope antibody (MCA 1360 from AbD 
Serotec). The activities of RlioA, Racl or Cdc42 were measured 
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Figure 7. Analysis of bone length of the Arhgap2Sf^^' mouse. 
Length of tibia and femur measured from X-ray images of 10 week-old 
wild type (n = 5), heterozygous (n = 9) and Arhgap28-null (n = 5) mice. 
doi:1 0.1 371/journal.pone.01 07036.g007 
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24 hours after transfection using G-LISA Activation Assays 
(Cytoskeleton), according to manufacturer's protocol. 

3D tissue constructs and recoil assays 

Primary embryonic chick fibroblasts were isolated from E14 
chick metatarsal tendons. For the formation of 3D tissue constructs 
in vitro, cells were seeded into fibrin gels as described previously 
[11]. To modulate Rho signaling in vitro tissue constructs were 
washed with phosphate buffered saline (PBS) three times and then 
equilibrated for 30 minutes with serum-free medium and then 
treated for 30 minutes with final concentration of 5 nM lysopho- 
sphatidic acid (LPA; Sigma) in 0.1% (v/v) DMSO or 10 ftM 
Y27632 (ROCK inhibitor; Sigma) in 0.1'^, (v/v) DMSO, diluted 
in serum-free medium, all at 37°C, 5% CO2 in a humidified 
environment. For the recoil assay, the tissue constructs were 
treated for 2 hours and was then cut from one suture to allow 
contraction to occur. The constructs were imaged using a digital 
single lens reflex camera at a fixed focal point for 30 minutes at 10- 
second intervals. The length of the constructs in each image was 
measured using ImageJ software. The length was converted to a 
percentage of the original and the means for each experimental 
group were calculated. The mean values were then fitted to a 3- 
parameter exponential decay function using SigmaPlot (Systat 
Software Inc.). One-way ANOVA and a Dunnett's test were used 
to determine significance differences between the derived T1/2 and 
recoil baseline values compared to control constructs incubated in 
DMSO. 

RNA and PCRs 

RNA was isolated with TRIzol reagent (Invitrogen) and treated 
with DNase (Promega) following manufacturers' protocols. cDNA 
was synthesized using TaqMan Reverse Transcriptase reagents 
(Applied Biosystems) and used for analyses by PGR and qPCR. A 
list of primers used can be found in Table S2. All PGR products 
produced from primers were validated and confirmed by DNA 
sequencing using the BLAST software against the NGBI 
nucleotide database and aligned with the expected sequence. For 
qPGR analyses, the 2"'^'^'^' method [48] was used to analyze 
relative fold changes in gene expression compared to the 
experimental control group. A two-sample t-test or one-way 
ANOVA followed by a Dunnett's test were performed to 
determine significance differences compared to the control sample. 

Immunofluorescence 

Gells fixed with freshly prepared 4'% (w/v) paraformaldehyde in 
PBS for 20 minutes at RT, then washed with PBS containing 0. 1 % 
(v/v) Tween 20 three times and blocked with 1% (w/v) bovine 
serum albumin, 0.1% (v/v) Triton X-100 in PBS containing 0.1% 
Tween 20 for 1 hour. Cells were incubated with primary 
antibodies, mouse monoclonal anti-V5 epitope or the appropriate 
control IgGs, diluted in the blocking buffer overnight at 4°G. Cells 
were washed with PBS containing 0.1% (v/v) Tween 20 for 10 
minutes three times and incubated with Alexa Fluor 594- 
conjugated antibodies (Invitrogen) and Atto 488-conjugated 
phalloidin (Sigma) diluted with the blocking buffer for 1 hour at 
room temperature, protected from light. Stained cells were washed 
with PBS containing 0.1% (v/v) Tween 20 for 10 minutes three 
times and mounted using Vector Shield containing DAPI (V ector 
Laboratories). Fluorescent images were taken using a digital 
camera attached to an Olympus BX51 and captured using 
MetaVue imaging software (Molecular Devices). 



Reporter and knockout mice 

Using homologous recombination in agouti C57BL/6N em- 
bryonic stem (ES) cell, the LlL2_Bact_P targeting gene trap 
cassette (including the genes encoding |3-galactosidase and 
neomycin) was introduced into intron 6 of the Arhgap28 gene 
by the NIH Knockout Mouse Project (KOMP; CA, USA). In 
addition, exons 7-9 are flanked by loxP sites. ES cells containing 
the Arhgap28 gene trap allele {Arhgap28^) were injected into 
C57BL/6J blastocysts to produce germ line-transmitting chimeras 
(prepared at the University of Oulu, Finland). Chimeric males 
were assessed on coat color and mated with wild type C57BL/6J 
females to produce heterozygous progeny, which were then mated 
to ohtain Arhgap28^ mice. For the generation ofArhgap28 deleted 
exons 7-9 (Arhgap28^') mutants, male Arhgap28^ mice were 
mated with females from a deleter Cre transgenic mice (a gift from 
M. Briggs, University' of Newcasde, UK) to ablate the loxP-flanked 
exons 7-9 in all tissues. Fl offspring heterozygous for the knockout 
allele were then mated to generate knockout {Arhgap28''"') mice. 

Genotyping 

Mice were genotyped using genomic DNA extracted from ear 
punches of adult mice or amniotic sacs of embryos using 200 Hg/ 
ml Proteinase K (Invitrogen) in buffer containing 17.6 mM N- 
lauroyl sarcosine, 100 mM NaCl and 5"/o (w/\') Chelex 100 resin 
(Bio Rad Laboratories). Genotypes were determined using specific 
primer pairs (see Table S2) using BioMix Red PGR reagents 
(BioKne). The annealing temperature for the wild type/mutant 
allele PGR was 53°C. The annealing temperature for the 
Arhgap28''''' allele was touch-down from 70-60°G for 10 cycles 
followed by 20 cycles at 60°G. The annealing temperature for the 
Cre transgene was 50°C. 

Whole mount X-gal stain 

For detection of beta-galactosidase expression, E10.5-E15.5 
embryos were fixed for 1 hour in 3.7% (w/v) PEA in PBS pH 8 at 
room temperature. For older embryos the skin was removed 
before being fixed for 1 hour in 3.7% (w/v) PEA in PBS pH 7.4 at 
room temperature. Embryos were then washed with PBS 
containing 0.1% (v/v) Triton X-100 for 15 minutes twice. Each 
embryos was then stained with 20 ml of freshly prepared X-gal 
staining solution (1 mM X-gal (Qiagen), 5 mM potassium 
ferricyanide (K3Fe(CN)6), 5 mM potassium ferrocyanide 
(KiFe(CN)6), 1 mM MgClj in PBS containing 0.1 "/» (v/v) Triton 
X-100) for 24 hours at 37°C. After staining, the embryos were 
rinsed with PBS and post-fixed in 4% (w/v) PFA in PBS overnight 
at 4°G. Aft(;r post-frx, the embryos were either processed for 
paraffin embedding or dehydrated in 70% (v/v) ethanol for 
6 hours and cleared with glycerol (10 ml of each 30, 50 and 80% 
(v/v) glycerol in 1% (w/v) potassium hydroxide and then 100% 
glycerol), incubating at 37°C for 2-3 days each. 

Histology 

Sagittal sections (6 ^m thick) were cut from wax embedded 
embryos and mounted for haematoxylin and eosin staining. Slides 
were stained using an automated stainer and cleared into Histo- 
Clear (Thermo Fisher Scientific). For Alizarin red staining, 
sections were then washed with distilled water and stained with 
2% (w/v) Alizarin red pH 4.2 for 10 minutes, washed with distilled 
water and dehydrated. Images were captured using a Carl Zeiss 
Axiocam Colour CCD camera with associated AxioVision 
software. 
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Comparative gene expression microarrays 

For the comparison of bone tissues between wild type and 
Arhgap28''"^ mice tibia and iibula w(-re dissc-cted from three PO 
neonatal litters each from defin(;d brixxiings. Bones were 
incubated in 1000 U/ml bacterial coUagenase type 4 (Worthing- 
ton Biochemical Corporation) in 0.25% (w/v) trypsin (Invitrogen) 
for 25 minutes at 37°C with agitation every 10 minutes. The bones 
were then removed of any excess muscle and cartilage tissues and 
washed in PBS. RNA was isolated using a dismembrator as 
described previously [49]. Integrity and measurement of total 
RNA was performed using Agilent 2100 Bioanalyzer (Agilent 
Technologies). RNA was amplified by two-cycle cDNA synthesis, 
and then labelled cRNA was synthesized and hybridized to Mouse 
Genome 430 2.0 GeneChip arrays (Affymc-trix). Microarray data 
are available in the ArrayExpr(;ss database (www.ebi.ac.uk/ 
arrayexpress) under accession number E-MTAB-2296. Micro- 
array data sets were analyzed by dChip (DNA-Chip) Analyzer to 
normalize the array readouts [50]. Normalized readouts were 
analyzed using the Robust Multichip Average method as described 
by [51]. Principal component analysis (PGA) was employed to 
confirm that different variables were present as a rjuality control 
for the arrays, p values for each probe set were generated by 
Limma t-test and q values were subsequently generated by 
applying false discovery rate correction. Gene ontology analysis 
was performed on probe sets that have detected fold changes 
greater than 2 using Database for Annotation, Visualization and 
Integrated Discovery (DAVID) online tool [52]. 

Bone length measurements 

X-rays of mice were produced using a Flaxitron x-ray specimen 

radiograph system (Flaxitron Bioptics) and x-ray film (GE 
Healthcare). Bone measurements were taken from scanned 
radiographic images using ImageJ soft\\arc'. One-way ANOVA 
followed by a Dunnett's test were used to determine significance 
differences compared to wild type. 

Supporting Information 

Figure SI Cloning of a V5-tagged Arhgap28 over expres- 
sion construct. A. Diffc-rcncc in C-terminal ends of the two 
major Arhgap28 isoforms from amino acid residue 679. PGR 
strategy to determine the sequence of endogenous Arhgap28 
transcript expressed in primary mouse fibroblasts. Overlapping 
products of PGR reactions 1 to 12 by specific primers (numbers in 
grey) were sequenced. Reactions in the boxes are unique to either 
the larger or the smaller isoform. ATG and TAG indicates the 
start and stop codons, respectively. Diagram not drawn to scale. B. 
RT-PGR products with RNA was isolated from primary mouse 
fibroblasts. A cDNA clone of the smaller Arhgap28 isoform was 
used as control template. C. Arhgap28 was cloned into a pcDNA6 
vector where a V5-His6 tag was introduced into the Arhgap28 
sequence immediately 5' of tiic stop codon. D. RNA was isolated 
from NIH3T3 fibroblasts transiently transfected with Lipofecta- 
mine only (L), empty vector (E) or Arhgap28-V5 clones (#6 or 
#8). RT-PGRs were performed to detect the expression of house- 
keeping gene, GAPDH (221 bp) or Arhgap28-V5 using specific 
primers (329 bp). No RT controls confirmed the absence of 
plasmid DNA contamination. E. Protein was isolated for western 
blotting to detect V5-tagged protein (■xpr(-ssi()n. A positive signal 
was detected at ~90 kDa. The predicted molecular weight of 
Arhgap28-V5 is 85 kDa. Gontrol lysates were included - Talin-V5 
(60 kDa), mock transfection (— ) and myristoylated-FAK-V5 (-I-; 
175 kDa). 
(TIFF) 



Figure S2 Arhgap2S^ mice are normal and express 
Arhgap28 due to unsuccessful gene trapping. A. Schematic 

showing the genot)'ping strategy for identifying the presence of the 
gene trap cassette targeted to the Arhgap28 gene. B. DNA was 
isolated from wild type, Arhgap28^'^' and Arhgap28^"'^' mutant 
neonatal tail tendons cells to confirm genotypes. DNA from wild 
type (+/+) animals wiU only produce a 493 bp product whereas 
DNA from homozygous (gt/gt) animals will only produce a 354 bp 
product and DNA from heterozygous {+/gt) animals will produce 
both bands. C. RNA was also isolated and RT-PGR was 
performed to detect the expression of Arhgap28 transcript 
spanning from exon 6 to 11, Collal and Gapdh. D. Sagittal 
sections of wild type {+/+), heterozygous {+/gt) and homozygous 
Arhgap28'^'mutant igt/gt) embryos at gestation day E15.5 stained 
with H&E and close-up of the limbs. 
(TIFF) 

Figure S3 Arhgap28^''' mice express an Arhgap28 tran- 
script lacking exons 7 to 9. A. Schematic showing the 
genotyping strategy for identifying the presence of the gene trap 
cassette targeted to the Arhgap28 gene and for the detection of 
Arhgap28 del7-9 KO allele. B. Representative gel image of 
genotyping PGR products. In the first genotyping PGR which 
distinguishes between wild type (493 bp) or mutant Arhgap28 
allele (either the Arhgap28^' or Arhgap2^'^ allele; 354 bp). The 
second genotyping PGR tests for the presence of the Arhgap28''''' 
allele, the product of the mutant Arhgap28 allele after Cre 
recombinase-mediated DNA excision (400 bp). The third geno- 
typing PGR tests for the presence of the deleter Cre transgene 
(350 bp). C. RNA was isolated from wild type, Arhgap28'^''''' and 
Arhgap28f^"'''^"^ pups and RT-PGR was used to detect expression 
of wild type Arhgap28 (634 bp) and Arhgap28^'' (338 bp) 
transcripts spanning from exons 6 to 11. RT-PGRs for Collal 
and Gapdh was used as loading controls. D. Sagittal sections of 
wild type and homozygous Arhgap28^^ embryos at gestation day 
E15.5 stained with H&E. 
(TIFF) 

Figure S4 Translated sequences of Arhgap28 tran- 
scripts. Amino acid sequences transcripts of Arhgap28 from 
wild type, predicted Arhgap28^"\ actual Arhgap28^ bone tissue. 
Sequence in blue is a result of remainder sequence from Gre- 
mediated recombination. Sequence in green is the RhoGAP 
domain. 
(TIFF) 

Figure S5 Quality control of microarray comparing 
wild type and Arhgap28^''' bone tissues. A. Agarose gel 

showing total RNA isolated from triplicate samples of PO tibia and 
fibula tissues from wild type and Arhgap28^'"^^ mice. B. PGA 
mapping of variability between the array samples. C. dGhip 
analysis of microarray data sets. (1) Triplicate samples of RNA 
from tibia and fibula of PO Arhgap28^'^ (wt) and Arhgap28^''"'^'^ 
(del) mice indicated by the number at the end of each array name. 
(2) Median intensity of microarray chip for each triplicate of each 
experimental group. (3) Present (P) call percentage indicates the 
percentage of total probe sets detected. (4) Array oudier 
percentage indicates the percentage of probe sets that have 
oudiers in the average readout profile of each probe within a probe 
set. (5) Percentage single outlier indicates the percentage of probes 
that do not have the same intensity pattern of other probes within 
the probe set. D. Gene expression changes between PO wild type 
and Arhgap28^'^ bones. Two-dimensional hierarchical cluster heat 
map analysis of the microarray readouts for all probe sets. This 
type of clustering is based on similarities in the expression profiles 
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and expression levels. Hierarchical clustering was performed using 
Cluster 3.0 software and visualised using Java Tree View. 
(TIFF) 

Table SI Summary of studies in which Arhgap28 is identified as 

a candidate of interest. 

(DOCX) 

Table S2 PGR primers. 
(DOCX) 
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